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Abstract. As a part of the dynamical heating emphasized by Hines, the contribution of 
aerodynamically produced atmospheric acoustic waves to the warming of the polar night 
mesosphere is re-examined with the polar night jet stream in the stratosphere as the source, 
and with the thermal structure of the polar upper atmosphere playing a part in propagation 
and absorption. It is shown that owing to the relatively small transmissivity of the thermo- 
sphere for long-period sound waves, most acoustic heating does not take place above the 
200-km level, but rather below 100 km, or around the mesopause, and that although acoustic 
output power from the polar jet stream is more than 10 s times larger in winter than in sum- 
mer, acoustic heating is not sufficient to compensate for the cooling rate around the meso- 
pause in winter, unless the wind velocity of the polar night jet stream continually exceeds 
2°° m/sec. 





Introduction 

An effective heating mechanism for the polar 
mesosphere, which is warmer in winter than in 
summer, is Kellogg’s chemical heating through 
the recombination of atomic oxygen [. Kellogg , 
1961; Young and Epstein , 1962; Maeda, 1963]. 
Although the velocity of subsidence required 
for this process is significantly less than that 
required for adiabatic heating alone, it must 
still be of the order of 0.1 cm/sec near the 60- 
to 90-km level. 

On the other hand, theoretical calculations by 
Haurwitz [1961] indicate an upward rather 
than a downward flow in this region. And al- 
though R. Sawada (private communication, 
1961) has questioned Haurwitz’s omission of 
certain inertia terms in his basic equations, a 
recent analysis by Kochanski [1963] of rocket 
data relating to mesospheric circulation shows 
good agreement with Haurwitz’s results. 

Rocket soundings by Nordberg and Smith 
[1963] have occasionally measured very warm 
wintertime mesosphere temperatures even in 
the middle latitudes where subsidence is un- 
likely, but where very strong westerlies in the 
region 40-80 km appear to be associated. 

Hines [1962] has suggested that the increase 
of temperature in the polar night mesosphere 
could be explained by dynamical heating. Gos- 
sard [1962] and Blamont and De Jager [1961] 
have found evidence of upward propagation of 
tropospheric disturbances into the ionosphere. 


The possibility of acoustic heating of the 
ionosphere was discussed by Daniels [1952], 
who considered noise produced by the sea sur- 
face as an energy source. However, Eckart 
[1953] showed that pressure waves from this 
source are completely ineffective in transferring 
energy up to such high altitudes. This has re- 
cently been confirmed by Cook and Young 
[1963]. The amount of acoustic heating of the 
ionosphere was also estimated by Golitsyn 
[19616], using tropospheric sources. He arrived 
at an upward energy flux of the order of 10"* 
erg cm"* sec, which is insufficient to cause the 
observed heating. Moreover, he did not consider 
the reflection of acoustic energy due to the 
vertical temperature distribution in the atmos- 
phere. \ 

In this paper, the contribution of aerody- 
namically produced atmospheric acoustic waves 
to the warming of the polar night mesosphere 
is re-examined, with the polar night jet stream 
in the stratosphere as the source, and with the 
thermal structure of polar upper atmosphere 
playing a part in propagation and absorption. 

It should be noted that the largest part of the 
kinetic energy of the atmosphere is in the tropo- 
sphere, where the air density is greatest, and 
is of the order of 10 s ergs cm If this tropo- 
spheric kinetic energy is not trapped effectively 
within the troposphere, the earth’s atmosphere 
at ionospheric heights should reach very high 
temperatures of the order of 10 s °K. However, 
Charney and Drazin [1961] have shown that 
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the escape of large amounts of planetary wave 
energy from the tropospheric layer into the 
upper atmosphere is prevented throughout most 
of the year by the large westerly zonal wind 
systems above the tropopause, except for short 
periods in the spring and autumn. 

Kinetic Energy Density of the 
Polar Statosphere 

Owing to the polar night vortex and jet 
stream, the kinetic energy of the polar strato- 
sphere increases remarkably during polar nights. 
As may be in the case of the solar corona, some 
part of this kinetic energy may be transferred 
upward by means of internal pressure waves. 
According to Boville et al. [1961], the kinetic 
energy density of zonal wind systems between 
60 and 80°N at 25 and 500 mb increases sig- 
nificantly during the wintertime (from Decem- 
ber until March). The kinetic energy density 
corresponding to the winter peak of the basic 
zonal flow is of the order of 150 and 600 ergs 



Fig. 1. Monthly average of kinetic energy den- 
sity at 30 mb (solid lines) and at 50 mb (dashed 
lines) during the period from July 1957 to July 
1958 above Ft. Churchill (59° N geographic lati- 
tude) in ergs cm'*. 


cm"* at the 25- and 500-mb levels, respectively. 
For example, the increase of kinetic energy at 
Ft. Churchill during the period of July 1957 to 
July 1958 (Figure 1) [U. S. Weather Bureau , 
1961] in the upper stratosphere in wintertime is 
caused mainly by the development of the polar 
night jet stream [. Krishnamurti , 1959]. 

The dissipation of the energy of jet streams, 
as well as any strong wind system in the atmos- 
phere, is mainly due to friction in the boundary 
layers, which cause turbulent motion to decay 
into small eddies in the surrounding mediums. 
These eddies decay further into thermal motion 
and cause local heating. However, some part 
of the kinetic energy in turbulent motion in jet 
flow can propagate into the upper atmosphere 
in the form of pressure waves. Owing to the 
exponential decrease of air density with height, 
the amplitudes of those propagating pressure 
waves grow rapidly in high altitudes and ap- 
proach shock-wave magnitudes; i.e., dissipation 
becomes nonlinear. 

The propagation and absorption properties 
and energies of these internal atmospheric waves 
are discussed mathematically in the next section. 

Equations of Internal Atmospheric Waves 

The notations and constants . 

a 2 , thermal conductivity of air, i.e., a 2 = k/C p p 
in cm 2 / sec. 

c, sound velocity in the air, km/sec. 

D/Dt, the Eulerian derivative, i.e., d/d^+V-V. 

f, resultant of all external forces. 

< 7 , acceleration of gravity, 9.8- 10 -3 km seer 2 

Hj scale height of isothermal atmosphere, i.e., 
H = RTo/g. 

Jo, Joule's mechanical equivalent of heat, 
4.185* 10 7 ergs/cal. 

k, horizontal wave number, corresponding to 
X, i.e., k = 27 r/X, in km -1 . 

p, p, T y small departure from static value of 
pressure, density, and temperature, function 
of x and z. 

Poy Po, To, static pressure, density, and absolute 
temperature which are the functions of z 
only in dynes/cm 2 , g/cm s , and °K, respec- 
tively. 

p, p, Ty total pressure, density, and absolute 
temperature, i.e., p = p 0 + p, p = p 0 + p, 
and T = T 0 + T. 

R, gas constant of air, 2.87* 10“ 4 km 2 sec~ 2 
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Fig. 2. Propagation surfaces (heavy lines) of internal atmospheric wave in the isothermal 
atmosphere with scale height 8 km ( T 0 = 0°C), given by (12); heavy lines and dashed lines 
correspond to sonic waves (acoustic waves) and therm obaric waves (internal gravity waves), 
respectively. 

V(u, w), velocity vector, where u is horizontal 
(southward) and w is vertical (upward) com- 
ponent of air motion in m/sec. 

z, period of pressure wave. 

7 , ratio of specific heats of air C p /C 9 , where 
C p and C, are the specific heat of air at 
constant pressure and that at constant 
volume, respectively. 

k, heat conductivity of air in cal/cm sec °K. 

X, horizontal wavelength of pressure wave, in 
km. 

fiy viscosity of air, in dynes sec/cm 5 . 

v , kinematic viscosity of air, i.e., v = fi/po 
in cm 2 /sec. 

r, period of pressure wave in sec. 

x(x, z), the divergence of velocity in sec -1 , 
i.e., x = du/dx + dw/dz where x and z are 
horizontal (southward) and vertical (up- 
ward) coordinate, respectively. 

co, angular frequency of pressure wave corre- 
sponding to r, 2t/ r. 

The fundamental equation. The equation of 
internal atmospheric waves is derived from the 
following three fundamental equations: namely, 
the equation of motion 

DV 1 

^ + -Vp + g = f (l) 


the equation of continuity 

£(i)- 1 , vv (2) 

and the equation of thermodynamics 

Dp/Dt = c(Dp/Dt ) (3) 

w 



Fig. 3. Diagnostic diagram of internal atmos- 
pheric waves in the isothermal atmosphere with 
scale height 8 km (To = 0°C) ; ordinate and ab- 
scissa are angular frequency w in sec” 1 and hori- 
zontal wave number k in km -1 , respectively. 
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Fig. 4. Diagnostic diagram of internal atmospheric waves in the isothermal atmosphere 
with scale height 8 km (T 0 — 0°C), where ordinate abscissa are period r in seconds and 
horizontal wavelength X in kilometers, respectively. r A and tb correspond to the wi and wb 
of Figure 3, respectively. 


Since the source of disturbance is extended 
along a certain latitude, the problem can be 
treated in two dimensions and the coordinates 
specified within a meridional plane. Let x be 
horizontal (southward) and z vertical (upward). 
In the first approximation, the Coriolis force and 
external forces, except gravity, are neglected. 
The equation of motion 1 then becomes 

Po(du/dt) — — dp/dx 

and 


p Q (dw/dt) = —dp/dz — gp 

The equation of continuity 2 for the first-order 
approximation is 

dp | dp 0 / .\ 

_ +w _ = _ Po . x (4) 


elimination of u, w, p, and p: 



Since c 2 = ygH, this is written further as 


* + -!<//-- nf* + 

dzH dz 


Lt gH 


- k 


. JiJL 
+ c C 2 H 



where H f = dH/dz. 

The diagnostic diagram. The wave equation 5 
can be written 


where x = x( x t z > 0 is the velocity divergence 
defined in the preceding section. By using (4), 
we can write the equation of thermodynamics 3 as 

( dp/dt ) + (c 2 x — wg)p 0 = 0 

The waves traveling in the flat atmosphere 
can be found by assuming that u, w, p } and p 
are proportional to exp (c ct — kx ). Then, by 


^ - 2N & + M\ = 0 (6) 

dz dz 

where 

N = (l/2fl)(l - H') (7) 

M 1 = ^ - k 2 + % (h' + 1 - (8) 

c a ) H \ yj 
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Fig. 5. Atmospheric resonant sound period ta 
and Brunt-Vaisala period tb as a function of alti- 
tude for the polar atmosphere in summer. 



PERIODS (sec) 


Fig. 6. Atmospheric resonant sound period ta 
and Brunt-Vaisala period tb as a function of alti- 
tude for the polar atmosphere in winter. 


Solution of the differential equation 6 is given 
approximately by 

x(x, z) = e N ‘(Ae~ n ‘ + 2fe + "') (9) 

where A and B are constant and can be deter- 
mined by boundary conditions. 

It should be noted that the solution 9 of (6) 
or (5) can be classified as cellular and non- 
cellular, 1 depending on the relative magnitudes 
of N and M ; i.e., 

1. The cellular solution: 

n = iv rj 2 = M 2 — N 2 for M 2 > N 2 

( 10 ) 

2. The noneellular solution: 

n 2 = N 2 - M 2 for N 2 < M 2 (11) 

The cellular waves have propagation vectors in 
both the horizontal (k) and vertical (rf) direc- 
tions, while noneellular waves have no propaga- 
tion vector in the vertical direction and simply 
propagate in the horizontal direction. 

Since M is a function of the horizontal wave 
number k f (10) gives the relation between two 


1 The nomenclature of cellular and noneellular 
is taken from Martyn’s [19501 paper, and these 
terms correspond to internal and external waves 
discussed by Chamey and Drazin [1961]. 


wave numbers in the vertical (rf) and in hori- 
zontal direction ( k ) : 


2 2 
rjc 


CO 

where 




+ 


7 2 2 / 2 2 \ 

k C {(j) 0) B ) 

2 / 2 2 \ 

CO (co — C C A ) 


( 12 ) 


or 


or 



- 55 <> - 

■’-Msr + <*-»») 


(14) 



co a and 00 B ar> e called the atmospheric resonant 
sound angular frequency and Brunt-Vaisala 
angular frequency, respectively [ Tolstoy , 1963]. 
These frequencies are both functions of the at- 
mospheric temperature distribution in the verti- 
cal direction and are related to the stability of 
layers in the atmosphere. 

We can see from (12) that the propagation 
surface is an ellipsoid if co > coa, but it is a 
hyperboloid if co < co#. In the former case the 
waves are called sonic or acoustic waves and in 
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Fig. 7a. Square of the index of refraction K 2 for atmospheric acoustic waves with period r 
in the summertime polar atmosphere as a function of height, where the parameter X is the 
horizontal wavelength in kilometers and r = 100 sec. 


the latter thermobaric [E chart, 1960] or internal 
gravity waves [Hines, 1960] (Figure 2). In both 
cases, the rotation axis is vertical ( z axis) . 

Since rj 2 becomes negative for o B < to < co^, 
no cellular wave exists for this frequency range 
(see Figures 3 and 4). These figures are called 
diagnostic diagrams of atmospheric waves 
[Eckart, 1960], in which a domain indicated by 
n 2 > 0 corresponds to the noncellular mode and 
two-domain n 2 < 0 (i.e., rp > 0) stands for the 
cellular modes. 

Approximate Transmissivity of the 
Polar Atmosphere 

Since the noncellular (or external) waves, 
whose periods are between t a and r B} have no 
vertical propagation vector, these waves neither 
contribute to vertical energy transfer nor exist 
locally in free atmosphere. For cellular (or 
internal) waves, on the other hand, the ther- 
mobaric (or internal gravity) mode can propagate 


with a period larger than Brunt period t b , but 
sonic (or acoustic) modes exist only below a 
period t a . 

As can be seen from (13) and (14), these two 
critical periods are a function of the atmospheric 
temperature and the vertical temperature gra- 
dient in the atmosphere (i.e., lapse rate of 
atmospheric temperature). Therefore, these 
values must differ not only with height, but also 
with season (Figures 5 and 6). 

The atmospheric models for the two seasons 
are essentially the same as the ones reported by 
Stroud et al . [1960] [see also Maeda, 1962, Figure 
13; Young and Epstein , 1962, Figure 1], In a 
domain between two critical periods, i.e., r A 
(solid line) and r B (dashed line), only the non- 
cellular mode can exist. Therefore no ther- 
mobaric waves with a period shorter than 370 
sec can propagate across the 50-km level and no 
sonic waves with a period longer than 250 sec 
can propagate through the 80-km layer in 
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Fig. 7b. Square of the index of refraction X* for atmospheric acoustic waves with period r 
in the summertime polar atmosphere as a function of height, where the parameter X is the 
horizontal wavelength in kilometers and r = 10 sec. 


summer. In winter both waves are able to 
propagate through these layers and the critical 
periods approach each other: 320 sec for ther- 
mobaric w r aves and 290 sec for sonic waves. 

To see the transparency (or transmissivity) 
of the atmosphere for internal waves more 
quantitatively, the wave equation 6 should be 
transformed into a more simple form. By sub- 
stituting for z 

y = f dz/H 
J *0 

and 

X® = e' /2 <Kv) 

the term dx/dz in the wave equation can be 
eliminated. Then (6) is written as 

<?<t>(y)/dy 2 = — K*<t>(y) (15) 

where 

X* - H 2 


As can be seen from (7), (8), and (10), tj s is a 
function of w, k , H , and H * . 

The one -dimensional wave equation 15 indi- 
cates that, if K 1 > 0, the solution co(y) has a wave 
nature, and corresponding waves can propagate 
along the vertical direction, within the domain, 
where this condition is satisfied. On the other 
hand, from Figure 2 we can see that there is a 
shortest (critical) horizontal wavelength \ e for a 
given period of the waves which corresponds to 
a largest horizontal wave number k c , and that 
the wave propagating close to vertical direction 
has longer horizontal wavelength than the one 
propagating in an oblique direction with the 
same period of the wave. In Figures 7 and 8, 
therefore, the waves propagating close to vertical 
axis can propagate upward as far as r < r Aj 
while the waves propagating obliquely meet a 
barrier, where, even if the condition r < t a is 
satisfied, K * becomes negative. If there are more 
than two barriers, those obliquely propagating 
waves are trapped between these barriers, form- 
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Fig. 8a. Curves similar to Figure 7 a for wintertime polar atmosphere; r = 100 sec. 


ing a duct in the horizontal direction. By 
comparing Figures 7 and 8, we find the barriers 
to be weaker in winter than in summer. 

Using the diagram shown in Figure 2, for the 
conditions at each level the critical angle 6 C can 
be drawn as a function of altitude for a given 
period r of the waves, provided that r < r A . 
6 e is defined such that the waves propagating 
obliquely with zenith angle larger than B c cannot 
enter into the above layer, simply because K 2 
becomes negative beyond this angle. In other 
words, 6 e can be regarded as the aperture of 
allowed cone around the vertical axis at each 
level in the atmosphere. Therefore the relative 
transmissivity / t (r, z) of each layer for vertical 
propagation of atmospheric acoustic wave is 
given approximately by(20 c /7r)*, provided that 
6 C is expressed in radians. From Figure 9, it 
is seen that the atmosphere above the meso- 
pause, which is nearly 80 km in the polar region, 
is always opaque for the atmospheric acoustic 
waves. It should be noted that the wind system 
in the higher layers is not affected very much 


by the relative transmissivity, because an oblique 
axis tilted windward can be regarded as the 
direction of maximum propagation. 

Upward Energy Transfer by Pressure 
Waves 

Acoustic waves generated aerodynamically in the 
atmosphere. Sound waves generated aerodynam- 
ically are different from sound waves produced 
by the vibration of solids. A strong jet flow in 
the free atmosphere produces sound waves due 
to turbulent motions in the boundary layers 
between the core of a jet and surrounding air. 
Lighthill [1951, 1954] investigated this problem 
and derived a theoretical expression for the 
power of sound waves as a function of jet stream 
velocity. His theory was further extended by 
Proudman [1952] to the generation of sonic 
waves by isotropic turbulence. 

Although the acoustic heating of the solar 
corona is discussed by Schatzman [1961], similar 
calculations cannot be applied to the earth's 
atmosphere, because the situation in the earth's 
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Fig. 86. Curves similar to Figure 76 for wintertime polar atmosphere; r = 10 sec. 


upper atmosphere is entirely different from that 
of solar atmosphere. As was discussed by Kidsrud 
[1954], the effect of magnetic fields on the 
generation of sound noise by isotropic turbu- 
lence is important as well as on the dissipation 
processed in the case of solar corona. 

Since the aerodynamical production of acoustic 
waves has been developed in detail by Lighthill 
and Proudman, here only their results will be 
applied to the present problem in the earth’s 
upper atmosphere. The acoustic power output 
per unit mass of air P, in ergs/g sec, is given by 

P = aeM 5 for M « 1 (16) 

where a is a nondimensional constant, c is the 
mean rate of energy dissipation per unit mass, 
and the Mach number M — u/c (u and c are 
the mean velocity of jet flow and sound velocity, 
respectively). 

According to Batchelor [1956, p. 103], e is 
given by 

e(3= du/dt) = Au/l 


where A is a nondimensional constant of the 
order of unity, l is a characteristic length, which 
can be taken as the reciprocal of the wave 
number at which the maximum of the energy 
spectrum occurs, and u is the root mean square 
velocity of turbulence. 

If the double correlation function given by 
Heisenberg [1948] is correct, the constant a in 
(16) is numerically 38. Then (16) can be written 

P ££ (70u 3 /Z)(M) 5 

In the upper atmosphere, where the Reynolds 
number is very large, the spectrum of turbulence 
consists generally of three ranges: 

1. The energy-producing range, which con- 
tains large eddies. 

2. The inertial subrange in which energy is 
neither created nor destroyed. 

3. The energy-dissipating range, which con- 
tains small eddies. 

As Proudman [1952] pointed out, the turbu- 
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Fig. 9. Relative transmissivity of polar atmos- 
phere, jt(r, z), for atmospheric acoustic wave with 
a period r = 100 sec; solid line and dashed line 
correspond to wintertime and summertime, respec- 
tively. 

lence which makes an appreciable contribution 
to the generation of sonic waves at large Rey- 
nolds number belongs to the eddies that do not 
dissipate energy. In other words, the main source 
of aerodynamically produced sonic waves is 
turbulent motion around jet flow, which belongs 
to (1) the energy maintaining range, and to 
(2) the inertial subrange. Therefore, according 
to Batchelor [1956, chapter 7], the following 
assumption can be made: 

10 3 < l < 10* and 1/10 < u/u < 1/3 

The distribution of turbulence along the jet 
stream can be assumed uniform ; then the acoustic 
output per unit area along the jet J(u), in 
ergs/ cm 2 sec, will be 

7.10"VM s po < l(u) < 3.10~VM 5 Po (17) 

where p 0 is mean density of air around the source. 

From Figure 10 the production of sonic waves 


is seen to decrease drastically below a wind 
velocity u of the order of 50 m/sec. Since the 
derivation of (17) is based on the assumption 
that the Mach number is small [ Lighihill , 1954], 
calculations with the output power beyond 
u = 150 m/sec might be incorrect. 

Attenuation of atmospheric waves. As has been 
shown, the propagation of atmospheric sonic 
waves is nearly isotropic if the wavelength or 
period is very short, while thermobaric waves 
propagate almost horizontally. However, if the 
wavelength (or period) increases, transmissivity 
in the vertical direction decreases. On the other 
hand, any atmospheric waves are attenuated in 
air by the following methods: (1) friction due to 
viscosity of air, (2) thermal conduction of air, 
(3) radiation through the atmosphere, (4) dis- 
persion due to the inhomogeneity of air in the 
atmosphere. Neglecting effects of (3) radiation 



Fig. 10. Acoustic output power of the polar 
stratosphere J(w), in ergs/cm 2 sec, with an as- 
sumed source height of 25 km, plotted against the 
mean wind velocity of the jet stream u in m/sec. 
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Fig. 11. Attenuation factor /„(r, z ) of atmospheric acoustic wave versus altitude z for 
t = 1, 10, 30, and 100 sec. Solid lines are calculated by the classical formula 18 and dashed 
lines are given by Golitsyn [1961a], taking Joule loss in the ionosphere into account. 


and (4) dispersion, which are generally small for 
acoustic waves as compared with the effects of 
(1) friction and (2) thermal conduction, we get 
the so-called classical Stokes-Kirchhoff formula 
for the attenuation coefficient a(r) [Schrodinger, 
1917; Rayleigh , 1945], which is given by 

«(r) = % a 2 + (18) 

rc i y 6 J 

where a* = K/C v p and v = p/p. The heat 
conductivity of air k is of the order of 6 X 10”* 
cal/ cm sec °K at NTP, and the specific heat C v 
is nearly 0.27 cal/g °K below the mesopause and 
increases with height to 0.3 cal/g °K above 
approximately 350 km [Ishikawa, 1959]. The 
viscosity of air p is of the order of 1.7 X 10“ 4 
poise at NTP. By using these numerical values, 
we can write the attenuation coefficient (equa- 
tion 18), in cm- 1 , 

a(r) ^ (2.4 X 10" 1 8 )/pr 3 


As Golitsyn [1961a] showed, long-period pres- 
sure waves undergo an additional dissipation 
owing to Joule loss of the current induced by 
fluctuations of the conducting medium of the 
ionosphere in the earth’s magnetic field. 

In Figure 11, the attenuation factor 

/„(r, z) = exp — £ a(r, z') dz'J (19) 

is shown for four different wave periods. 

An estimate of the rate of energy dissipation 
per unit time per unit thickness per unit cross 
section in the upper atmosphere is made by the 
following approximation: 

= f z) -j / 0 (w, z) dto (20) 

where dEjdz is in ergs/cm 3 sec, and / a (co, 2 ), 
the attenuation factor, is given by (19) and 
shown in Figure 11. /,(oj, z) is relative transmis- 
sivity of atmosphere as shown in Figure 9. The 
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Fig. 12. Relative dissipation rate of energy flux transferred by atmospheric acoustic waves, 
( dEt/dz)/Io in %/km, where I 0 — Io(u) is the acoustic output at source as shown in Fig- 
ure 10. Solid lines and dashed lines correspond to wintertime and summertime polar atmos- 
phere, respectively; heavy and light lines stand for flat spectrum {m — 2) and steep spectrum 
(m — 5) of the waves, respectively. 


differential frequency spectrum of the wave i(oi) 
will be assumed to be 

i( co) = iooT” 1 m > 1 (21) 

From Figures 5 and 6, the lower limit of the 
integral in (20), o) A , is to be taken as follows: 

Summer: 


CO A = 

w A (z) 

for 

0 < z < 20 km 


0.023 

for 

20 < z < 70 km 



for 

70 < z < 85 km 


0.025 

for 

z > 85 km 

inter: 




0) A = 


for 

0 < 2 < 30 km 


0.024 

for 

2 > 30 km 


Strictly speaking, the right-hand side of (20) 
should be multiplied by a geometrical factor 
which is inversely proportional to the distance 
from the source. The width of the source, the 
polar night jet stream, is more than several 
hundred kilometers [Krishnamurti, 1959] and 
the boundary is widely diffused horizontally. 
For this reason, the geometric attenuation factor 
is not important and can be neglected in the 
present approximation. 

A constant to in (21) is given by 

/ co m — 1 

i 0 o)~ m do) = i Q — m > 1 

)A(* o> m — 1 

where 7 0 = Io(u) is the output power of acoustic 
waves at the level of source z 0 given by (17) 
and shown in Figure 10. 

In Figure 12, the upper two curves of the 
relative dissipation rate of energy (< dE t /dz)/1 0 
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Fig. 13. The rate of atmospheric heating dT(z)/dt, ° C/sec, in each layer due to the dissi- 
pation of acoustic energy flux h (ergs/cm* sec) versus altitude z. Solid, dashed, heavy, and 
light lines correspond to those in Figure 12. 


versus altitude are calculated by neglecting the 
transmissivity of the atmosphere; i.e., /,(«, z) is 
assumed unity. This corresponds to the esti- 
mation made by Golitsyn [19616]. 

If all those dissipated energies of acoustic 
pressure waves are assumed to be converted into 
thermal energies of air, the rate of heating at 
each level is given by 

dT(z) ^ 1 1 dE t 

dt C p p(z) Jo dz 

where p(z) is mean atmospheric density at the 
altitude z. 

The results shown in Figure 13 correspond to 
Figure 12 for a unit output power of the source, 


and, since the latter is shown in Figure 10, the 
following points can be made : 

1. If the wind velocity around the 25-km 
level is of the order of 30 m/sec, as in the usual 
summertime polar stratosphere, then, from 
Figure 10, the acoustic power at the source is 
less than the order of 10 -5 erg/ cm* sec. There- 
fore, the rate of maximum heating, dT/dt is of 
the order of 10“ 8 X 10“ 6 °C/sec or 10“ 8 °C/day 
at the maximum, of which the height is around 
75 km if the wave spectrum is flat (m = 2), or 
100 km if the spectrum is steep (m = 5), respec- 
tively. 

2. On the other hand, if the wind velocity 
around 25 km exceeds 100 m/sec (which is 
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steadily observed in the polar night jet stream 
[Krishnamurti, 1959; Riehl, 1962]), the output 
power h is between 10 3 4 5 and 1 er1 g/cm 2 see, as 
can be seen from Figure 10. The corresponding 
maximum warming rate is then 2 X (10 -5 ~ 10 -8 ) 
°C/sec or 2 X (1 ^ 10~ 3 ) °C/day at the altitude 
of 75 km for m = 2, or around 100 km for m — 5, 
respectively. 

3. If the transmissivity of the atmosphere is 
disregarded, as for Golitsyn’s estimate, then, 
from the two upper curves shown in Figure 13, 
the rate of maximum heating is of the order of 
lO-io °C/sec (10 -6 °C/day) in summer, and 10 -2 
to 10 -6 (10* — 1 °C/day) in winter. The height 
of maximum heating is around 260 to 300 km, 
corresponding to m = 2 to 5. 

Conclusion 

Although the available estimates of the 
acoustic energy flux for upper air heating have a 
wide range (Figure 10) due to ambiguities of 
several parameters, such as the mean velocity 
of turbulent motion around the jet stream, the 
effective size of eddy in this turbulence, and the 
frequency spectrum of turbulently produced 
acoustic waves, the following remarks can be 
made from present calculations: 

1. Owing to the temperature gradient, the 
atmosphere above the 100-km level is not trans- 
parent for atmospheric acoustic waves, especially 
long waves. Therefore acoustic heating is most 
effective for warming the 100-km level rather 
than 200- and 300-km levels (Figure 13). 

2. On the other hand, if the transmissivity of 
sound waves due to thermal gradient in the 
atmosphere is disregarded, most heating occurs 
around the 300-km level, and its rate exceeds 
10°C/day in winter. 

3. Although the acoustic power from polar 
night jet stream is more than 2 orders of magni- 
tude larger than that of a tropospheric disturb- 
ance, acoustic heating around the mesopause is 
not sufficient to compensate for the cooling rate 
(about 10°C/day) of these levels in winter. 

4. In this respect, another mechanism of 
dynamical heating such as energy transfer by 
means of internal gravity waves [Hines, 1963] 
might be effective for such warmings of the 
mesosphere as those observed by Nordberg and 
Smith [1963]. 

5. If the wind velocity of the jet stream 
continually exceeds more than 200 m/sec, the 


layer around the mesopause can be warmed 
significantly even by acoustic waves. 

Note added in proof: The relative transmissiv- 

ity of atmospheric acoustic waves, ft (r, z ), is 
more accurately given by 

/,(r, z) = 1 — cos d e (r, z) 

where 


and k c (r, z ) is the maximum horizontal wave 
number at the height z for a given period r, and 
z 0 is the height of minimum atmospheric tempera- 
ture. 

However, there is no essential difference in the 
general trend of ft (r, z) with altitude z as shown 
in Figure 9, even if more rigorous transmissivity 
is required. 
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